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The first successful macroscopic theory for the motion of superfluid he- 
lium was that of Lev Landau (1941) in which the fluid is modelled phe- 
nomenologically as an interpenetrating mixture of a superfluid and a 
normal fluid. It has later been shown that Landau's two-fluid model 
can be rigorously derived from a one-fluid model within the classical 
fields approximation. Assuming a separation of scales exists between the 
slowly varying, large-scale, background (condensate) field, and the short 
rapidly evolving excitations, a full description of the kinetics between 
the condensate and the thermal cloud can be obtained. The kinetics 
describes three-wave and four-wave interactions that resemble the C12 
and C22 terms, respectively, in the collision integral of the ZNG theory. 
The scale separation assumption precludes analysis of the healing layer 
and thus does not include the dynamics of quantised vortices. Whilst 
the analysis required the use of small parameters arising from the scale 
separation assumption and the assumption of a weakly depleted conden- 
sate, we expect the results to hold true over a wider range of parameters. 
This belief is motivated by the validity of Landau's two-fluid model which 
can be derived from a one-fluid model using nothing more than Galilean 
invariance principles. Indeed, we argue that similar arguments can be 
used to recover a two-fluid model directly from a classical field simply by 
invoking a local gauge transformation. This derivation does not require 
any small parameters to be introduced suggesting that the results that 
lead to the kinetic equations may turn out to be more general. 
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1.1. Introduction 

Historically, the development of the two-fluid model of superfluidity for liq- 
uid He H has relied on a phenomenological approach. By starting with the 
governing equations of motion for a single fluid, and introducing an addi- 
tional velocity field associated with the superfluid component, Landau [TJ[2] 
formulated his seminal two-fluid model of liquid He H using little more than 
the principle of Galilean invariance. Landau's model predated the discovery 
of quantized vortices which were conjectured more than a decade later by 
Onsager [3] and Feynman [J. These conjectures received further support 
with the theories of Gross [F and Pitaevskii |6] formulated for a weakly 
interacting Bose gas that were derived from microscopic considerations. 

Some authors, notably Hills and Roberts J^, Ginzburg and Pitaevskii 
[8], Khalatnikov [9 , and Geurst [10], have sought generalizations of Lan- 
dau's two-fluid model by including relaxation effects as well as the effects 
of healing. The latter is particularly important in light of the presence 
of quantized vortices. In fact, it is well established that the Euler equa- 
tions, governing the motion of an ideal classical fluid, preserve the topology 
imprinted by the initial conditions onto the flow. In contrast, vortices in 
superfluids undergo reconnections (e.g. ^llj ) which consequently have im- 
portant effects on their dynamics [12] . These can only occur if the effects of 
healing are retained in the model. The rationale for including the effects of 
healing into Landau's two-fluid model is, therefore, perfectly understand- 
able. However, the validity of such a two-fluid description for Helium II 
per se has been brought into question (Ginzburg and Sobyanin [E]). This 
is due to the healing length in He II being of the order of the inter- atomic 
distances thereby precluding a continuum description on these scales. 

To account for the effects of healing and relaxation in a physically consis- 
tent manner, we will review a more fundamental derivation of the two- fluid 
model first proposed by Putterman and Roberts [111 [15] (hereafter PR). 
PR arrived at the two-fluid model through an intermediate step involving 
the derivation of a kinetic set of equations. Their work is in some sense 
the basis of the classical fields approach. This was developed further by 
Kagan and Svistunov [TB], and Davis et al. [T7|, as a model for studying 
finite-temperature effects in atomic Bose-Einstein condensates. Our pri- 
mary goal in this chapter is, therefore, to review some of these results and 
to establish their connection to other closely related theories. In particular, 
we will show that the local-gauge transformation used by Coste [I^ is a 
limiting case of the PR theory. 
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1.2. Kinetics in Nonlinear Classical Fields 

In recent years, the classical fields approximation has evolved from a qual- 
itative method to describe the kinetics of condensate formation from a 
strongly nonequilibrium initial condition (BerlofF and Svistunov |19|). to a 
method that is used to model finite temperature Bose gases [2^. In the 
context of modelling superfluid He II, PR were two of the earliest advocates 
of the classical fields approximation. The key objective of these authors at 
the time was to elucidate the kinetics that ultimately give rise to a two-fluid 
description of liquid Helium II. Remarkably, by starting from the equation 
of a single fluid, and using a scale-separation argument, they were able to 
recover a set of kinetic equations governing the thermal excitations. In the 
collision-dominated regime this leads to the two-fluid model of Landau. 

Following Putterman and Roberts [21], we will use the GP equation in 
the classical fields approximation as our starting point. This will facili- 
tate in relating the results presented here to the application of local gauge 
transformations. The GP equation for the classical field ■ip{x,t) is 



If we define ip{x,t) = \^/j{x,t)\exp[iip{x,t)m/h], where p{x,t) = rn|?/'P, and 
u(x, t) = V</3, this equation can be written in hydrodynamic form as 



At this point, we wish to stress that the hydrodynamic form given by Eqs. 
p.2p must not be attributed to the superfiuid density and the superfiuid 
velocity. All we have done is simply to re-express the equation for the 
classical field in terms of an amplitude and a phase. 

To proceed, we will consider small amplitude excitations on top of a 
background field. We do this by assuming that the system contains two 
disparate length and time scales. We can then decompose the field vari- 
ables into a long-wavelength slowly varying background field denoted with 
subscript '0', and short wavelength rapidly varying excitations, denoted 
by primed quantities. Formally, we introduce a slow timescale r and a 
long-wavelength scale X such that t = t6 and X = x5. Then we have 



<p(x,t) = (^o(X, r) +.^'(X, r,x,t), p(x, t) po(X, r) + p'(X, r, x, t). (1.3) 



The small amplitude requirement is enforced by requiring \ip' /ipol ^ 0{e) <C 
1 and likewise for \p' / pq\. Moreover, given the scale separation assumption. 



Zm 



(1.1) 




(1.2) 
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we will decompose the perturbation into 

(p' = e[(^i(X,r)+J^2(X,T)]e^«, p' = e [pi(X, r) + (5p2(X, r)] e*®. (1.4) 

The perturbation is, therefore, represented by a monochromatic wave with 
a varying wavelength controlled by the parameter S ^1. With d0(x, i) = 
k • dx — u)dt corresponding to the rapidly varying scales, we have, 

a;(X, r) = ~dQ/dt, k(X, r) = VO. (1.5) 

We shall begin by assuming that <C (5 <C e <C 1. We then seek 
how excitations propagate on some initial background state denoted by 
/9(X,0) = /9o(X) and v(X,0) = Vo(X). Substituting Eqs. ((O]), into 
Eq. (|1.2I) and expanding to 0(e), we obtain the local dispersion relationship 

a;=Zj + k-vo, ZU± = ±co|k| = ±cofc, cq ^ iy/gpo)/m. (1.6) 

where uJ is the intrinsic frequency, (the frequency in the frame of reference 
of Vo) and Cq is the speed of sound. At ©(eJ), we obtain a conservation law 

dA d 



-{Cg,A)^0, (1.7) 



for the waveaction per unit volume defined hy A = S/TD, £ = Cglpip/po = 
Pok'^lfi]'^ , where Cgi = doj/dki is the group velocity. This conservation 
law has been obtained by considering a monochromatic, slowly varying 
wavetrain, of characteristic wavenumber k = k(X, r). A general initial 
perturbation, pi(x, 0) say, can be represented by a linear superposition of 
such wavetrains which evolves into pi(x,t) at time t. In view of the scale 
separation that is assumed to exist in our system, we introduce the Gabor 
transformed quantities, such that 

/oo 
p'(x',t)e*«-(''-''')^(|x-x'|)dx', (1.8) 
-oo 

with T^(X,k, t) defined in an analogous way from ip' . Following [2?, we 
take = {S* /V2^)^ exp{-{S*^)^ /2} and set 5 < ^* < 1 so that the 
kernel T varies on the intermediate scale. The inverse transform is 

pi(x,0) = ^-^=^ j i3(ko;X,0)e'''°-^dko, k(ko; X, 0) = ko, (1.9) 

pi(x,^) ^ / B(k;X,r)e^Q('';-^*)dko, k(ko; X, r) = VO. (1.10) 

(V27rd*)-^ J-oo 

The mean waveaction per unit volume can now be defined as 

^(k; X, r) = (co/pofc) {\B{k; X, t)\') , (1.11) 
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with the waveaction density per unit volume given by the integral of 
yl(k;X,T) with respect to the initial wavenumbers ko. To express this in- 
tegral over the wavenumbers k, we introduce n = aA where a denotes the 
Jacobian 9ko/9k. This Jacobian satisfies Liouville's theorem (Soward p3] ) 

do- da dcgi 

Combining the above Eqn. with Eqn. (|1.7I) and then using the chain rule 
to transform from the independent variables (ko,X, r) that is used in Eq. 
(|1.7p to the variables (k, X, r), we finally arrive at the wave Vlasov equation 

dn dn du dn /-, 

As is clear from the form of the equation, this governs the motion of the ex- 
citations in the collisionless regime. We note that n(k, x, t) can be identified 
with Bogoliubov quasiparticles (see e.g. Lvov et al. [22]). 

Now in the opposite regime, where we have 5 ^ 5* ^ ^ 1, we 
obtain, at leading order, a kinetic equation with a collision integral arising 
from the nonlinear terms. Under these conditions, the kinetic equation can 
be derived in a rigorous way using the methods of matched asymptotics. 
Since these derivations are quite involved, we will merely sketch out the 
initial key steps and defer the details of the closures used to the relevant 
references. In what follows, we will suppress the dependence on X since 
it does not play a key role in the kinetics in the parameter regime we are 
considering here. It then follows from Eq. (|1.2I) and the above definitions 
that to 0{e^) we have 

^+ik-vo^ B{k) - pk^V{k) = e Jji-k2B{ki)V{k2)dki. (1.14) 

Similarly, we can derive an equation governing V"(k). By introducing 

B{k, t)=J2 Oe"'""'*, ^(<9a7at)e-*'^-* = 0, (1.15) 

we can derive from Eq. (jl.l4p and the corresponding equation for V{k) an 
equation for the evolution of the Fourier amplitudes a* (k, t) given by 

= E / ^(ki,k2)a^Hki)a^Hk2)e^(^^-^^^^-^=^^)*rfki. (1.16) 

Si ,S2 ki+k2=k 

The coefficient H is defined in PR. Starting from the above equation, we can 
derive equations for the second order correlation which would correspond 
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to the waveaction density n(k). However, in so doing, we find that it 
depends on the fourth order correlation that arises in the integrand. The 
equation for this fourth order correlation in turn depends on even higher 
order correlations and so the equations cannot be closed. Benney and 
Saffman [21] (see also derivations by Connaughton and Pomeau PS|) show 
how the hierarchy can be closed under an appropriate set of assumptions 
(e.g. quasi-Gaussian statistics). Further details of the derivation can be 
found in PR but it can be shown that we finally obtain the kinetic equation 



If we extend the above analysis to the next order in e, then wc find 



X [n(ko) + n(ki)] - n(ko)n(ki)[n(k2) + ^(kg)]} 
X(5(w2 +UJ3-UJ0- ZUi)(5(k2 + k3 - ko - ki)dkidk2. 

The coefficient K is defined in PR. The first term proportional to in Eq. 
()1.17|) represents resonant triad interactions of the waves with wavenumbers 
ko,ki,k2. The resonance conditions contained in the delta functions arise 
from conservation of momentum and conservation of energy of the gov- 
erning equations. We note that this leading order behaviour of the wave 
kinetics, derived for a condensate containing a large fraction of the total 
number of particles, has a similar structure to the C12 term in the collision 
integral of the Zaremba, Nikuni, and Griffin (ZNG) model j26]. On the 
other hand, the four- wave resonant interactions, described by Eq. (jl.lSp . 
corresponds to the C22 term in the ZNG theory and represents interactions 
of quasiparticles within the the thermal cloud. In fact, the key difference 
between the classical field and ZNG model stems from the missing sponta- 
neous scattering contributions that are retained in the ZNG theory. Given 
that we have started from a mean classical field representation which as- 
sumes macroscopically occupied modes (f^(k) 3> 1), it is natural that a 
classical field cannot model spontaneous scattering processes. 

Using the above results, we can now recover a hydrodynamic two-fiuid 
model. To help identify the microscopic basis of the results to be derived 




(1.17) 



-n(ko)[n(ki) + Ti{k2)]} 6{ujo - ^1 - uj2)S{ko - ki - k2) 
2|H(-ki,k2)P {n(ko)n(ki) ~ n(k2)[n(ko) + n(ki)]} 





(1.18) 
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in the next section, using local gauge transformations, we have opted to 
present a brief derivation. We recall that in the collisionless regime, corre- 
sponding to <C (5 ^ e <C 1, we obtained the Vlasov wave equation. On the 
other hand, the kinetic equation was obtained in the collision-dominated 
regime when 5 <C e'* <C 1. If we relax our condition on 6 by requiring 
e,6 <^ 1, then we expect to obtain a Boltzmann like equation of the form 

— +Cg-Vn-V^-Vkn=(—) . (1.19) 



dt " \dt J_ 

From this equation, we can show that 

d f f f duj 

— j nojdk -|- Vj / nujCgjdk— / n-^dk, (1-20) 

— J nkidk + Vj J nkiCgjdk = —J n-^—dk, (1-21) 

I](n)dk-FVj J ^{n)cgjdk>0, (1.22) 

where S(7i) is the entropy density. In the collision-dominated regime, we 
have local thermodynamic equilibrium with {dn/ dt)co\\ — 0. The equilib- 
rium distribution can then be written in terms of the temperature T as 

n = neq(/317), (3 = l/kBT, 17 = w-k-w. (1.23) 

We note that the equilibrium Rayleigh- Jeans distribution given above is the 
long-wavelength (classical field) limit of the Bose-Einstein distribution. We 
also note that n — ricq is the only solution for which equality is obtained 
in Eq. ()1.22|) . This follows from the H-theorem corresponding to the case 
where the entropy is maximised. Now using neq, we can define the following 
thermodynamic variables at each point X, 

En ^ J ricqUJdk, p„w = J Ticqkdk, S = J Y,{ncq)dk. (1.24) 

Hence, from Eq. (|1.21l) and (ll.22|) . we obtain 

d d dT dv - 

— + V • (5v„) = 0, (1.26) 
at 

where the entropy per unit volume [S) is advected by the normal fluid 
component. To obtain the equation of continuity and superfluid velocity, 
we expand the original equations of motion up to 0{e'^) so that 

^ + V • (povo + p„w) = 0, ^ + vo • Vvo = - Vm, (1.27) 



January 28, 2011 1:7 



World Scientific Review Volume - 9in x 6in 



chapter2 



8 



H. Salman, N.G. Berloff, and P.H. Roberts 



where fj, is the chemical potential and is related to the pressure P by [3T] 



Attributing po;Vo,w to the total density p, the superfiuid velocity and 
the relative normal/superfluid velocity (v„ — Vg), respectively, we recover 
Landau's two-fluid model. We note that one can include dissipative effects 
by extending the above analysis using a Chapman-Enskog [27 ^ procedure. 

We end this section by asking why two-fluid phenomena, whose premise 
is a one-fluid model, cannot be observed in classical fluids. It turns out 
that the reason for this is set by the requirement p„ < p. This condition is 
typically satisfied only by He II at very low temperatures thus preventing 
this phenomena from being observed in other classical fluids. 

1.3. Two-Fluid Model from Local Gauge Transformations 

Akin to Landau's original derivation of the two-fluid model from a one- 
fluid model, PR have revealed the kinetic basis of the two-fluid theory 
and how it emerges from a one-fluid model. In that sense, their results 
generalise Landau's result since a kinetic description remains valid even 
when a hydrodynamic description breaks down. This occurs, for example, 
very close to T = where the mean free path of the excitations becomes 
so large that a hydrodynamic description of the normal fluid component no 
longer holds. However, one can also arrive at the two-fluid model directly 
from the GP equation. To show this, we recall that the GP equation can 
be derived from an action principle with the Lagrangian density given by 



The above Lagrangian density is invariant under the global gauge transfor- 
mation — > i/jc*'-"'"/''-' (a is a real constant) and under spatial and time 
translations of the form x — ^ x-|-Jx and t — )■ t+St. It follows from Noether's 
theorem that these symmetries lead to the conservation laws of mass con- 
servation, momentum conservation, and energy conservation, respectively. 
Since two-fluid hydrodynamics is a manifestation of an additional macro- 
scopic degree of freedom, we expect the respective conservation law can 
be attributed to another symmetry in the Lagrangian density. We argue 
that the new macroscopic degree of freedom is associated with a broken 
local gauge symmetry. Early attempts to derive two-fluid hydrodynamics 
from local gauge transformations were carried out by Chela-Flores [28], 





(1.28) 




January 28, 2011 1:7 



World Scientific Review Volume - 9in x 6in 



chapter2 



Two-Fluid Model of Superfluidity 



9 



and Cummings et al. [29^. However, the resulting equations were inconsis- 
tent with Landau's model since they did not respect Galilean invariance. 
Coste [TB] presented an alternative derivation of Landau's two-fluid model 
using local gauge transformations. Coste proceeded in a phenomenological 
way resting his ideas on the assumption that the GP equation only describes 
the condensate close to T = 0. However, we argue that the approach does 
in fact have a more fundamental basis arising as a limiting case from the 
kinetic description presented in the previous section. In fact, in the limit 
where local thermodynamic equilibrium holds, the macroscopic variables, 
such as superfluid and normal fluid velocities, become functions of space 
and time. The local gauge transformation essentially involves introducing 
a local gauge field that is related to these variables as we will demonstrate. 

We begin by noting that the GP equation is invariant under a Galilean 
transformation of the form x — > x + Vi, and t — > t provided the wavefunc- 
tion is transformed according to V'(x, t) t) exp[i(— V-x-|-iF^<)m/fi]. 

Rewriting the field t) in terms of an amplitude jV'l and a phase (ipm/H), 
we have V(/3 — > 'Vip + V. The transformation, therefore, relates V' in one 
frame of reference to V' in a frame of reference moving with velocity V. If 
we demand that V(/3 correspond to the local superfluid velocity, then we 
must transform to a local frame of reference such that this requirement is 
satisfied. This is readily achieved by introducing a local gauge transforma- 
tion ip — >■ ^e*"^'''*^™/''. If we now define C = —dta and A = Va, we can 
determine the forms of ( and A by noting that for the Galilean transfor- 
mation given above, we have V — >■ V and dt — > 9t — V • V. Since the 
spatial derivative does not change under a Galilean boost, and A is defined 
as the spatial derivative of a, it follows that A must also remain invariant. 
Similarly, using the transformation of the time derivative with V coinciding 
with the local normal fluid velocity v„, it follows that 



X and T] are Galilean invariant scalars that are functions of the density and 
entropy only. Now under the local gauge transformation, the Lagrangian 
density given in Eq. (|I.29p transforms as 



A x{p, S){vs 



Vri), C ^ 1]{p,S) +Vn ■ A. 



(1.30) 



Ci{ip,i;*,S, v„) 



2,1 1 




Tn\i,\\ 



+ 



(1.31) 
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The Euler-Lagrange equation for now reads 

1^ + V • [p{l - xf^s + PX(1 - X)v„] = 0. (1.32) 
Equation (|1.32p becomes the two-fluid equation for mass conservation when 



Pn ^ px{p, S)[2 - x{p, S)], Ps = p[l - x{p, S)]"^ . 
Similarly, the Euler-Lagrange equation for p gives 



(1.33) 



P = V^ 



dt 2^ ^' 
drj 
Pd-p 



r 



9P_ 

m 2 



P- 



2m? 



(Vp) 

Ap 8p2 



21 



(1.34) 



These expressions show that x — > 0, and p — > ps, as T — > 0. It then fohows 
that the generahsed Lagrangian density reduces to the original GP equa- 
tion that governs the condensate motion close to T = when ( ^ p.. Thus 
far, we have recovered the equation for mass conservation, and the equation 
for the superfluid velocity. To obtain the remaining set of equations for the 
two-fluid model, we need to introduce additional constraints on the system 
through the use of Clebsch potentials in line with Geurst's [10] variational 
formulation. These additional constraints are required since it has been 
known for some time |30' 'ST that the equations of motion of hydrodynam- 
ics, when derived from a variational principle, tend to be constrained to an 
irrotational flow. This is associated with a particle relabelling symmetry 
that is lost in the Eulerian representation of a fluid and must be enforced as 
a separate constraint through the use of Lagrange multipliers. The distin- 
guishability of a 'fluid element' is associated with the thermal excitations 
that satisfy an independent equilibrium distribution at each point in space. 
The final form of the Lagrangian density is therefore given by 

C2=Ci+ a[dtS + V • (5v„)] + ^[dtiPS) + V • (/35v„)]. (1.35) 

with the two constraints corresponding to conservation of entropy and con- 
servation of relative fluid vorticity. Now on the hydrodynamic length scales, 
the quantum pressure term appearing in Eq. (jl.31l) can be neglected. In 
this limit we can recover Landau's two fluid model as described in 1181. 
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